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ON THE NATURE OF THE 
f RADIAL AND CROSS TRACK ERRORS 

4 

FOR ARTIFICIAL EARTH SATELLITES 

N. L. Bonavito 
R. A. Gordon 
J. G. Marsh 

ABSTRACT 

In this paper we discuss the analysis of the radial and cross track errors 
of artificial Earth satellites in terms of the interference of two one dimensional 
celestial mechanical wave trains. The resulting equations for these tracking 
errors describe the behavior of the uncertainties in the orbital parameters as 
oscillatory in nature, with a rapidly oscillating term, which is a function of the 
sum of the observed and computed orbital frequencies, modulated in amplitude 
by a slowly varying oscillation. This latter term is itself a function of either 
the difference between these orbital frequencies or between the values of the 
computed and observed right ascensions, depending upon whether it is the radial 
or cross track case under consideration. In theoretical physics, this effect is 
referred to as the phenomenon of beats and is common in the physics of sound 
and elastic media. These results indicate that the cross track calculation describes 
the behavior of uncertainties in the right ascension of the ascending node and 
the inclination, while th i radial calculation gives information on uncertainties in 

v 
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the semi-major axis, the eccentricity, and the argument of perigee. Analysis is 
done for the U. S, Air Force Cannonball (OAR-901) and GEOS-n satellites. 

In addition, expressions for the radial and cross track oscillatory fre- 
quencies are obtained in terms of the orbital frequencies of the satellites. These 
oscillatory frequencies are functions only of the zonal harmonic ♦ :ms of the 
Earth’s gravitational potential, and are used to analyze the behavior of the 
tracking errors in terms of uncertainties in the gravitational field coefficients. 

Finally, we show that the time average of the radial and cross track errors 
in any case, will both approach zero. 
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I. INTRODUCTION 

There has been considerable interest in the error analysis of spacecraft 
trajectory systems. Error analysis can be defined as the ability to describe 
the effect of inherent uncertainties of an overall orbital determination system 
on the computational accuracy of the system. There are several sources re- 
sponsible for the presence of such uncertainties, with perhaps the most im- 
portant due to the inability to correctly describe the physics of the problem. 

As a result, the physical data that is used, and the modelling of the forces involved, 
provide at best, only a good starting point for the calculation of an orbit. The 
accuracy of the observational data will depend upon the type and quality of the 
tracking technique that is employed. Even after fitting an orbit, the post con- 
vergence residuals do not account for uncertainties in coordinates of the track- 
ing sites. More important, even though the constants of integration may be well 
determined, the accuracy of the calculated orbit will still depend upon the accu- 
racy of the differential equations of motion. In any event, the overall uncertainty 
in the calculated coordinates of a spacecraft can be defined within the framework 
of a time dependent error bound (Reference 1). One method of studying the effect 
of uncertainties in orbit determination is to examine the radial and cross track 
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addition, the radial error is defined as this same difference projected along the 
instantaneous unit radius vector of one of the orbits. If the orbit system can 
be considered as a conservative one, then calculation of the radial and cross- 
track errors can be used to study errors produced by various gravity models. 
Such results are of considerable importance in the field of geodesy. 

In this paper we shall review these concepts in the light of classical theo- 
retical physics and ( ompare the results with real and simulated data of various 
orbit systems. 


n. STATEMENT OF THE PROBLEM 

Let us now define an inertial coordinate system, fixed at the Earth's center 
with the x-axis pointing toward Aries, z through the geographic North Pole, 
and y set so as to form a righthanded system. Then the observed position vector 
of a point mass in orbit is given by r 0 = T x 0 +7 y 0 + k z 0 , where T, J and k 
are the orthogonal unit vector set. The observed velocity is given by v 0 = i v x + 
jv y +kv i , and the computed values of these quantities are designated as r c = 
Tx c + Ty c +kz c ,andv c = "iv x + Jv + kv r respectively. From the 
above discussion we then have that the cross track error 

<» 

|E„I 

where L Q = (r Q x v 0 ), andAr = (r 0 - r c ). The radial error Is given by 

R = — • Ar (2) 
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Before proceeding with a theoretical analysis, let us briefly examine these two 
equations. If we assume that our system is a conservative one, that is, if there 
are no forces present drawing energy from the orbit, then we have that the 
Hamiltonian, namely, the kinetic plus potential energy, is the total energy. In 
addition, if time is not explicitly present, then the Hamiltonian is also a constant 
of the motion Furthermore, the total angular momentum and the z-component 
of angular momentum are constants of the motion. Now the gravitational potential 
energy can be written as a generalized solution of Laplace’s equation, namely, 


V 


tL 

r 



P*(cos 6) (C cos mX. + S sin m\) 

n ' ' n , m n , m 


where m is the product of the gravitational constant and the planet's mass, r e 
is its equatorial radius, r and 0 are the planetocentric distance and colatitude 
of a field point, C n m and S are the coefficients of potential, \ is the angle 
of east longitude, and P"* (cos 0) are the associated Legendre functions. In 
practice, one tries to include all forces in nature acting on the system in the 
Hamiltonian, and in the process of orbit determination, a complete knowledge 
of the geodetic constants and coefficients of the Earth’s gravitational potential 
is necessary for perfect orbit calculation. Since such information is not available, 
and in addition, since present computational techniques are insufficient to main- 
tain agreement with observation for more than a few days, then it is customary 
to compensate for these errors by reinitializing the constants of integration 
through a fitting procedure by comparing theory with observation over an arc of 
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several days. However, it appears that there may be a limitation to what we may 
know about the gravitational field. For example, the possibility of mass trans- 
port in the liquid outer core of the Earth, even if uniform, would still give rise 
to density fluctuations. In addition, mass transport arising from rapid first order 
phase transitions (condensation of matter under enormous pressure) within the 
lithosphere could create local density changes which would correspond to a 
creation and destruction of higher order coefficients in the gravity field. More 
advanced theories indicate that rapid first order polymorphic phase transitions 
may be the source of deep focus earthquakes, manifesting as irregularities or 
breaks in the Chandler wobble uf the Earth. It is also estimated that during 
such an event, enough mass is shifted so as to cause as much as a 0.1 to 0.4 
milligal change in the local value of gravity. In any event, such considerations 
seem to place a temporary if not permanent restriction on a complete theoretical 
description of the gravitational field and the subsequent accuracy of the pre- 
dicted positions of a spacecraft. 

In spite of such limitations however, we can still make the following 
rationalizations: Since the total energy, and the angular momentum are constants 
of the motion, then the energy difference of the computed and observed orbits 
which is manifest through the spacecraft i adius vector by way of the semi- 
major axis, is a bound quantity. As a consequence, one might expect the radial 
error to be bound, and since the angular momentum is conserved, the cross-track 
would be likewise. 


4 



/ 


7 

III. INTERFERENCE OF WAVES 

The problem that now suggests itself from the above discussion, is the 
nature of the tracking errors arising from an inability to properly or adequately 
specify the potential energy through the coefficients of the gravitational field. 

In actual practice, tracking errors are often calculated by one of the following 
methods: For a given potential energy, two partially overlapping arcs of fitted 
data are obtained on a spacecraft by comparing observations of position and 
velocity with those calculated using the given potential. The positions and 
velocities in each of the overlapping portions of the fitted arcs corresponding 
to a particular time are then used in equations (1) and (2) to determine the 
cross and radial track errors. In this case, the potential is fixed, and the 
coordinates (within the fitted arcs) are changed. In the other case, positions 
and velocities of two slightly different orbits are compared at corresponding 
times, and the tracking errors are again computed for each of these times. Here 
the initial conditions and the equations of motion are the same, but values of 
the coefficients of the Earth’s gravitational potential differ slightly for each 
orbit. 

Technically speaking, interference is a term referring to the physical 
effects of superimposing two or more wave trains. As in the case of elastic 
media, it is desirable to define a ’wave equation* which contains information on 
the behavior of the system. We now consider the polar equation of an ellipse, 
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r = a(l - e cos E) (3) 

where a, e, and E are the semi-major axis, eccentricity, and eccentric anomaly 
respectively. A one dimensional representation of this equation is obtained by 
multiplying by the cosine of the true anomaly. Therefore, 

y = r cos v = a(l - e cos E) cos v, (4) 

and using the anomaly connection for the true and eccentric anomaly, 

(cos E - e) /5\ 

cos v = ' , 

1 - e cos E 

y - a cos E - ae. (®) 

This is our one-dimensional wave equation for elliptic motion of an artificial 
satellite, referenced to the major axis of the ellipse. 

Let us now consider two waves of equal frequency and amplitude travelling 
with the same speed in the same direction, but with a phase difference <t> between 
them. Therefore, with the conditions that a. l = a 2 = a, ej ^ % , and E = 2^11, we 
have for the two waves 

Yj = a cos(27rf 1 t) - acj 
y 2 =: a cos(2rrf j t - £) - ae^ 


(?) 

( 8 ) 
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The resultant wave which is the difference of equations (7) and (8) is. 


y = y, - y 2 = (ac 2 - aCj) - a [cos (27j ^t) - cos(2^f % ♦ - ;:)] . 


From the trigonometric equation for the difference of the cosines of two angles. 


cos A - cos B = - I s i n i (A + B) • s i n i (A ~ B). 


we obtain, 


- a(e 2 - Oj) - pasin-jj sin 2 t; , t 


This resultant wave corresponds to a new wave having the same frequency f t 
but with an amplitude, 2a«sin^/2. If '£ is very small, the resultant amplitude 

m 

will be nearly zero. As a result, equation (10) suggests that the radial t acking 
error as given by equation (2) will oscillate with the same frequency as the 
orbital frequency but with amplitude determined by any phase difference between 
the computed and observed orbits. A phase difference in the eccemric anomaly 
however, is equivalent to one in both the true and the mean anomaly also. S ; nce 
the orbital periods are the same here, then a phase difference between the two 
eccentric anomalies arises for example from a poor determination cf the 
argument of perigee. Near the beginning of the orbit therefore, this uncertainty 
in the argument of perigee will be a factor in the magnitude ol the rapid 
oscillatory term of equation (10). 
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Let us now consider the case of two wavetrains similar to the above, but 
with slightly different frequencies, traveilirg through the same region. The two 
displacements can be represented as 

y x = a co;> 27/fjt - ae r (11) 


and 


The difference is 


y 2 = a cos(27rf 2 t - <£) - ac 2 . 


( 12 ) 




y = - V 2 = a(e 2 - Cj) - 2a si 


The resulting vibration may then be considered to have a frequency 


7 _ f . + f 2 


which is the average of the two waves, and an amplitude given by the expression 
in the brackets. Therefore, the amplitude varies with time vrith a frequency: 


5 = 


f I- f 2 


If f l and f 2 are nearly equal, this term is small and the amplitude fluctuates 
slowly. This phenomenon is a form of amplitude modulation and is commonly 
referred to as 'beats'. A beat, or a maximum of amplitude w«ll occur whenever 
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equals 1 or -1. In orbital motion, the frequency of rotation is governed primarily 
by the value of t* * semi-major axis. This demands that y* and y 2 of equations (11) 
and (12) contain different values of a, which would result in a numerical determi- 
nation of y t - y 2 versus time. However, since amplitudes only add together 
algebraically, the general behavior of the system is still described analytically 
by equation (13). The interference of two waves of different amplitudes, fre- 
quencies, a ?A phases, would ordinarily result in a distorted complex waveform, 
but the analysis of equations (1) and (2) can be made more directly through 
equations (10) and (13). Equation (13) implies that errors in the computed value 
of the semi-major axis, and to a lesser degree, the ^centricity and argument 
of perigee determine both the amplitude and the frequency of the envelope of the 
radial tracking error (Figure 1). Equation (10) indicates that any distortion in 
the envelope is due to a timing error in the mean anomaly. If the semi -major 
axes of the two systems under comparisons are almost equal initially, then the 
corresponding frequencies or periods are likewise almost identical. Because 
of this, the advance of perihelion for each orbit is approximately the same, and 
the angular difference between each perigee (or apogee) increases slowly. 


When this difference is tt, the respective apogees are exactly out of phase, and 
the magnitude of their difference is at its maximum. This is also the maximum 
of the amplitude of equation (13). If the periods are long, £u i in addition, if 
they are close in value, then the time between succeeding maxima may be very 
large. Figure 2 shows a plot of the radial error versus time tor the cannonball 
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| (OAR-901) satellite (a = 8335.34 km, e = 0.123, and i = 92.00 degrees). The 

; graph was generated using the Brouwer theory of artificial satellite motion 

t . • • 

(Reference 2). For the calculation, the potentials used only the first two zonal 
harmonics, with the second zonal coefficients given by J 2 = -1.0825 x 10~ 3 and 
J* = 0.90 x J a . in this case, the period between maxima (or minima) is approxi- 
mately 120 days (Figure 2). The value of the eccentricity determines the 
symmetry character of the oscillation. If e is much different from zero, then 
Kepier f s law demands that the spacecraft spend more time near apogee. When 
the two orbits are essentially coincidental, the por'^ionai differences are 
smallest As the angle between the major axes increases, these differences 
increase, and the maximum is obtained when this angle is v 9 and as described 
above, when the spacecraft are at their respective apogee points. Between 
7 t and 277, the position differences decrease, becoming a minimum at 277. If the 
reference orbit is ’larger* that is, if the semi-major axis of the reference 

; orbit is larger than that of the ’computed* orbit, then radial error will *arc* 

positively. As the angle between the positional differences, in the plane of the 
reference orbit increases, the magnitude of their difference will also increase, 
and the envelope of equation (13), given by 




will also reach its maximum value given by twice the semi-major axis. The 
time for this to occur will depend on the value of f x - f 2 . In Figure 2, the 
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maximum of the variation is approximately 1G.000 km. The constant term 
a( e 2 - e j) also represents a shift about which the radial error oscillates. 

We can now summarize briefly. A comparison of a plot of equation (2) 
and equation (13), shows the following: An 'error 1 or uncertainty in the 
eccentricity will appear as a shift along the ordinate axis, about which an 
oscillation in time of the radial error occurs. This oscillation is of frequency 
(f x +f a )/2 where f % and f 2 represent a computed and observed orbital frequency. 
This rapid oscillation which is close to both f x and f 2 is then modulated by a much 
lower frequency oscillation of value (fj - f 2 )/2, and with amplitude 2a. The 
frequency of the orbit varies inversely as the period, and the period is a function 
of the semi -major axis. As a result, a difference between the computed and 
observed values of the semi-major axes results in a finite value of f t - f 2 , and 
a slow modulation of the rapid radial oscillation appears. The result 
suggests a beat phenomena, common in the theory of sound. The beat frequency 
is given by (f x - f 2 )/2, and as fj approaches f 2 , that is, as the difference between 
the computed and observed values of the semi-major axes approaches zero, the 
period of the beat phenomena becomes large. 

We now turn to the discussion of the nature of the cross track error given 
by equation (1). Let us assume that at the initial time, the six orbital param- 
eters which describe the computed system differ slightly from those of the 
observed orbit. At some later time, we also assume that the component 
of the positior vector of the spacecraft from the observed orbit into 

11 
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the corresponding line of nodes is r 2 cos u 2 , where r 2 = a 2 (1 - e 2 cos E 2 ), 
u 2 - j 2 + is the argument of perigee, and u 2 is the true anomaly. For 

the computed orbit, we have similarly, r x cos \i l along its line of nodes. The 
wave trains, represented by r 2 cos u 2 and r 5 cos u t need to be referenced 
to a common axis to discuss interference effects. However, for the cross- 
track error component, we need the perpendicular component of the observation 
on the calculated orbit, and hence, we wish to have the two wave trains interfere 
with each other but at right angles. Physically, this is the same as sighting 
parallel to the r l cos u t wave, and observing the oscillatory behavior of the 
r 2 cos u 2 v'ave with respect to the first. Therefore, the projection of r 2 cos u 2 
onto the line of nodes of the computed orbit is given by r 2 cos u 2 -sin (0 2 - ), 

where 0 2 and are the respective right ascensions of the ascending node for 
each orbit. Since the observed system is Inclined at an angle i 2 , the projection 
is not perpendicular to the computet} orbit. The perpendicular component is 
obtained by a rotation through this inclination. The result is 

Yq = a 2 sin i 2 * s in(0 2 - Qj) [cos co 2 * cos E 2 

- /I - s in cj 2 • s in E 2 - e 2 cos . (14) 

However, since the observed and computed orbits are not parallel, there is 
another contribution to the cross-track term originating from the discrepancy 
in the inclinations. Results show that both this term and its variation are 
very small. Nevertheless, its contribution is. 
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Yi = r 2 sin(i 2 - i,) = a 2 (l - e 2 cos E 2 ) sin(i 2 - i,). (15) 

The sum of (14) and (15) give the result we want, nan o l y, 

y = yf/ + y» = a 2 s in i 2 • s in(0 2 - [cos cj 2 ■ cos E 2 
- vl - c 2 sin • sin E 2 - e 2 cos a> 2 ] 

+ a 2 (l - e 2 cos E 2 ) sin(i 2 - ij). ( 16 ) 

If e is sufficiently small, we have, 

y = a 2 [s in i 2 • s i n(0 2 - 0 1 )*cos(E 2 + cj 2 ) + s in(i 2 - ij)] . (17) 

The angular velocity of the eccentric anomaly is much greater than that of the argu- 
ment of perigee (usually two orders of magnitude), so that the cross-track error 
displays a rapid oscillation with a frequency approximately that of the orbital 
frequency, modulated by a slowly increasing sinusoidal function with maximum 
amplitude a 2 sin i 2 . This is again analogous to beat phenomena, and the 
respective maxima are obtained when the two orbit systems differ in their 
values of the ascending rode by v/2. Figure 3 shows a plot of the cross-track 
error for the Cannonball satellite using the above, values of J 2 and J*. The beat 
period or that period between successive minima (or maxima) is seen to be 
approximately 870 days for this case. Equations (16) and (17) suggest that as 
the values of the computed and observed right ascensions become more in 
agreement, the beat period becomes progressively larger, with a maximum 


13 



16 


amplitude of a sin i. In our example, i = 92°, sothatsini^r 1, and the maximum 
theoretical cross track amplitude according to equations (16) and (17) is approxi- 
mately 8000 km, which is in close agreement to that given by equation (1). 

Figures 4 and 5 illustrate the use of equations (1) and (2) in the determina- 
tion of the cross and radial tracking error for a small portion of the orbit of the 
GEOS-IE satellite (a = 7700.00 km, e = 0.032, and i = 105.80 degrees). The tracking 
data consisted of Navy Tranet Doppler observations recorded by a globally dis- 
tributed set of tracking stations. The orbits were computed using the GSFC 
GEM-1 (Reference 3) and the SAO 1069 (Reference 4) Standard Earth gravity 
models. Perturbations due to solar-lunar gravity, direct solar radiation pressure, 
and air drag were also modelled in the solutions, so that the primary error 
source is the Earth's gravity field. The orbital solutions were obtained using 
the GEODYN orbit and geodetic parameter estimation system. For both the 
GSFC GEM-1 and SAO 1969 cases, orbital arc lengths of two days were used 
covering the periods 12 hours, 15 minutes of May 23, 1968 to 12 hours, 15 minutes 
of May 25, 1968, and 12 hours, 15 minutes of May 24, 1968 to 12 hours, 15 minutes 
of May 26, 1968. The cross and radial tracking errors were determined during 
the one day overlap from 12 hours, 15 minutes of May 24, 1968 to 12 hours, 

15 minutes of May 25, 1968. Because of the short period over which the calcula- 
tions were performed, the periods of the beat phenomena for both the radial 
and cross track terms are not shown. According to equations (13) and (17), 
the corresponding maxima are 15,400 km and 7700 km respectively for these 


14 


17 


fitted arcs. An examination of the radial error for the GEM-1 case shows a 
slowly increasing envelope, but with a downward drift, indicating that the 
computed value of the semi-major axis is larger than the observed value. To 
. determine this latter value, we would proceed in the following manner: The 
radial oscillatory period, P is read off of the graph (Figure 4 or 5). From this, 


7 _ fp+ f c _l 

2 p 

Using the fitted data for this arc, the computed period f c is calculated by 


(18) 


s 

fc 


r = 2 v 


1/2 a 3/2 


(19) 


where a is obtained from the orbit improvement process- Therefore f c = 1/P C . 
The observed frequency f 0 is then given by 


and from this, 


f 0 = 2f - f e 


P - 1 
r o “ 

r o 


( 20 ) 


( 21 ) 


From the expression for the period such as equation (19), the observed value of 
the semi-major axis is obtained. The difference between this and the fitted 
value of the semi-major axis gives us the discrepancy in the latter. Figures 4 
and 5 are not sufficiently expanded to give accurate values of P and a (e 2 - ej). 
The magnitudes of the radial errors for these two graphs indicate a slightly 
more serious phasing proDlem or uncertainty in the argument of perigee for the 
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Standard Earth solution. On the other hand, the slope of the radial envelopes 
suggests that th^ Standard Earth value of the semi-major axis may be closer 
to the observed value than that of GEM-1. This could be resolved by a calcula- 
tion such as given above. 

For the cross track component, the time scale of figures 4 and 5 are too 
short to make an accurate estimate of the discrepancy in the right ascension of 
the ascending node for the Standard Earth and the GEM-1 cases, as described 
by equation (17). 

In this section, we have attempted to reduce the analysis of the cross and 
radial tracking errors given by equations (1) and (2) above, to an analysis of 
the resultant behavior of the interference of two o.ie -dimensional 'celestial 
mechanical wave trains’. These results, namely equations (13) and (17), then 
describe the behavior of, and the discrepancies between, the computed and 
observed values of the orbital parameters for a given potential or gravitational 
field for any arc of the orbit. The radial wave equation gives information on the 
semi-major axis, the eccentricity and to a limited extent, the argument of perigee. 

I ! 

r 

The cross track wave equation gives information on the right ascension of thr* 

i 

ascending node and the inclination. Thus, five of the six orbital parameters are 

! 

examined by this procedure. This leaves only the mean anomaly. To accurately 

j 

r assess its behavior and any discrepancy between the observed and computed value ; 

• j- 

of this parameter, it is necessary to make a careful analysis of the corresponding \ 
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in-track error component. From Kepler’s law, it is obvious that for different j 
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values of the sidereal mean motion, two orbits will sweep out di ferent equal 
areas and different values of the true anomaly in equal times. From the anomaly 
connections between the true and eccentric anomaly and using Kepler's law, the 
mean anomaly, 

M-E-esinE, (22) 

will vary accordingly. 

IV. FREQUENCY EQUATIONS 

In the previous section we discussed the analysis of the tracking errors by 
relating their behavior to uncertainties in the orbital parameters. By this 
approach, we wish to asceitain the degree to which these orbital parameters 
are involved in the accuracy of the trajectory calculation. In this section 

1 

however, we attempt to relate the behavior of the radial and cross track errors ; 

to certain of the coefficients of the Earth's gravitational field. As a result, some j 

I 

information can be derived from the relationship between these orbital parameters 

and the values of gravitational coefficients used in the calculations. \ 

As a starting point for our calculation, let us consider the following expres- j 

) 

sion for the gravitational potential V(r,6 y )i written in terms of an expansion in ) 

i 

spherical harmonics: i 

i 

f 
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In addition, we also consider the following problem: A spacecraft represented 
by a point mass moves in orbit about a planet whose gravitational potential is 
given by equation (23). At some time which we shall take to be the epoch of the 
orbit, an instantaneously small but finite impulse is applied to the point mass. 
The problem is to describe the subsequent motion. We note that the impulse 
has the effect of disturbing the coordinates of the satellite at this initial or 
epoch time, while the potential field is constant. The Hamiltonian for this 
system, which is also the total energy, can be written as 

& =s E = i mr 2 + i mr 2 # 2 + ~ mr sin 2 6<$ 2 + V(r, 6). (24) 

Here, the geocentric right ascension is cyclic, and as a result, its canonically 
conjugate momentum 

p^ •= mr 2 sin 2 &4> = constant. (25) 

We may introduce an effective potential f V f (r, 0) for the motiorn 

. . p| 

V (r. 0) r- V(r, 6) + (26) 

2mr 2 s i n 2 0 

If we now assume that the colatitude coordinate is held fixed, and there is a 
small change in r only, then we have from equation (24), 

E - 'V’(r) = ■- mr 2 . (27) 





/ 


-* v «*** T+r *•<*** i - 


21 


Since the right hand member cannot be negative, the motion must be confined to 
those values of r for which f V f (r) < E. The relationship between and the 
equilibrium radius r, , that is, the radius prior to the applied impulse is, 


p •"<«■? [ _ / dV(r) \ _ P* = 0 

L dr -L ' dr 'r? mrA sin 2 0 


£ c 

Combining equations (28) and (25) we have. 


(28) 


= Hi) 1 

r \ dr / 2 r 

\ L , mr,TSiir (/ 

*V c 


(29) 


From equations (24) and (28) the equilibrium energy is 


-T^P) 

r e 


+ V(r£> + 2 


(30) 


For an energy slightly larger than E,„ and an angular momentum p given by 
equation (28), the coordinate r will perform simple harmonic oscillations about 


the value r^ . If we set 


k - [d 2 , V , (r) | _ / d 2 V(r ) \ f 3 P» 

L dr 2 J \ dr 2 / mrjsin 


(31) 


8 


f 8 e 


Then for small values of (r - ), we can expand ’V'(r) in a Taylor series: 


*V(r) = E fe , + *-k r (r - r_,) 2 . 


(32) 
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The energy equation (27), now becomes 

E-E^Imf 2 + Ik (r-r..) 2 < 33 > 

This is the energy for a harmonic oscillator with energy (E - Er,) t coordinate 
(r - r,), mass m, and spring constant k r . Then for a small impulse, the radial 
component of the radius vc 3 tor will oscillate about the equilibrium value r^> with 
a frequency given by. 



It is of interest to compare and i 2 : From (29) and (34) we have, 


(34) 



In other words, \ e now have the frequency of oscillation of the perturbed radial 
component of the radius vector in terms of the undisturbed or equilibrium value 
of the orbital period as given by the motion of the geocentric right ascension, 
and the first and second derivatives of the gravitational potential energy. A 
similar argument for the colatitude angle 0 f yields the harmonic oscillator 
energy equation, 


E - Eg = 1 mr 2 <?' 2 * i k(<9 - 6 >-,) 2 , (36) 

with energy (E - ), coordinate (0 - 6U), mass mr 2 , and spring constant k 9 

given by 
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K = 


"d 2 'V* (vj 

'= 

9 , 

ft\ 2 V(0)'\ 

(rin 2 J.+ 3 cos 2 ') 

dv 2 

l dv 2 / 

!• sin; 1 cos.- ' \ dr* / 

i ( 


The motion of the geocentric right ascension in the undisturbed mode is given 
by, 


*2 -- 


l 


mr 2 s i r. d. , cos c* ( , 


r vH 


(38) 


while after the impulse is applied, the colatitude angular component of the radius 
vector oscillates about the equilibrium value c» , with a frequency given by 


> 

ilQ 

II 

1 ** 

1 * 

II 

1 /d^V (< < \ 

\ 

4» 

(s i n 2 2>v + 3 cos 2 f\ ) 

M i 

mr 2 

mA dr 2 1 

mr 2 

sin cos ca. 

\ A . 


- ' c 




The relationships between u, d 

“d»« 

is. 





A?- 


(40) 


Here, the frequency of the angular component is also given in terms of the 
equilibrium value of the orbital period, and the first two derivatives of the 
potential energy, but with respect to the colatitude. 

At this point, it is instructive to choose some representations for the gravi- 
tational potential given by (23) and calculate c* r and for each case. 
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Case I. Let us assume that the Earth is perfectly spherical, and homogeneous 
in density. The potential is then given by 


V = ii 


( 41 ) 


Here then, 


dV(r) ^ 



and 


d 2 V(r)_ 2 fj. 



and from equations (31) and (25), 


From (29), 


k r ^j = ItL + 3m sin 2 0<p 2 . 


which yields, 


and 


*r(M) 




mro sin 2 d 




_ 


CO 


2 

r (M) 


k. 

m 



(42) 


(43) 


(44) 


(45) 
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iD which case, 




“rOO = sin& * r(M) (4? 

Since dV(r)/d# = 0 for the central force potential, then co e (v) = 0. Equation 
(47) is an oscillation in r superposed upon a motion around the z-axis with an 
angular velocity given by equation (29); will vary slightly as r oscillates, . 

but will remain very nearly equal to the constant value given by equation (43)* 

If our reference frame is rotated so that the z-axis is always perpendicular to 
the orbital plane, then sin 0=1, and equation (47) reduces to 


r(M) V r(/1)‘ 


For the case of the spherical Earth, the radial component of the perturbed radius 
vector has the same period of oscillation or rotation as that of the orbit, and will 
execute one oscillation per revolution of orbit. 

Case II. Let us now consider the case where the Earth is represented as a 
perfect rotational ellipsoid, and its potential is given by 


V j^l + ~ (3 cos 2 6 - 1)J . 


Proceeding as in the above case we have for the rrdial component, 




i mm 


v\v:o- 



i 
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‘^r(M+J 2 ) 3 ~ 2 " [* + 2 (“l Jz(3 

2 mr , sin 2 t/ J * \ T pj 

C \ ' 


cos 2 Q - 1) 


( 50 ) 


and 


r<^+J 2 > 


m mr 2 
t 


1 "4 f-) Jt ( 3 cos 2 9-1) 


= sinfc i 


For the 0 component, 


m 


1/2 


(3 cos 2 e - 1) 


^£ 5^(3 cos 2 <9 - 


1) 


r < *W+j,r 


j, - 3Mja N? 


(51) 


(52) 


(53) 


, k a 3 M J 2 /r e 

^(m+J 2 ) = — = 3 

mr 2 mr J 


( t ) cos2 ® e - 


(54) 


and 


*&</*+ j 2 ) ~ cos ^0<^+j 2 r 


(55) 


Equations (52) and (55) are again referenced to our original fixed inertial z-axis. 
By rotating the coordinate system so as to make the z-axis lie perpendicular to 

t 

the orbital plane, we again have that sin B ~ 1, and cos 0 = 0, in equation (52), 


while a rotation of x-coincidental with the radius vector wiT make cos 0 o = 1 in 

c. 
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equations (55). A transformation of these quantities into the planes of motion of 


the tracking errors are necessary to describe the oscillatory behavior of the 
latter in terms of u> T and . The fact that = 0, for the spherical Earth 

results of course, from the fact that we have in that case only central Ijrce 
motion. As a result, if a small impulse were applied normal to the plane of 
motion, the effect would be o y to change slightly, the inclination of the orbit. 

It would however, remain fixed at the new inclination. There would be no oscillatory 
motion in d\ The second case, namely that of an Earth's gravitational potential 
containing a zeroth, and second zonal harmonic only is somewhat more interest- 
ing. Equation (55) where cos c? =1, indicates that the 0 -component executes 
one full cycle per revolution of the orbit, similar to the radial component given 
by equation (47). AfUr rotation, equation (52) reduces to, 


r(^ + J 2 > 


3 Ij 2 I rA 


3 iJ 2 !A, 


2 r, 




Here, we have included the sign of , the coefficient of the Earth's second zonal 
harmonic. Since the magnitude of J is approximately 0.10825 x 10~ 2 , then the 


coefficient of <£ r(/i + j > is only slightly less than unity. This means that the 
frequency of oscillation of the radial component of the perturbed radius vector 
is also slightly less than the rotational frequency of the spacecraft in its orbit. 
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From this then, the period of revolution of the radial component given by = 

r 

27 t/ca‘ t , will b i \u than the natural period of the orbit, =2 7 t/4' t . This 

frequency } uift, and correponding change in period is induced by the inclusion 
of the J 2 te ,t\ or second zonal harmonic in the Earth's potential. As a result, 
the mass redistribution is equivalent to a change in the spring constant of a 
translational mechanical circuit or the capacitance of the electrical series 
circuit. In this case, our results are now written, 




In words, if a small but finite impulse is applied to an Earth satellite in orbit 
about an Earth described by a perfect rotational ellipsoid, the radial component 
of the radius vector will oscillate about its equilibrium value with a frequency 


i 

whose value is i 




S5B?* ffSTTOVzr TT m ic » j w ui w w ip m m wm 1 r^gy-Mii rryrm f rr^rnr^n^y^ 


t 

l 

l 

f 

fc' 


ii 


times that of the geocentric right ascension. The frequency in 6 is identical to 

that of the geocentric right ascension. Also, the indication here is that we should 

expect a shift in both the frequency of the radial and cross-track errors, as 

either the number of, or the values of the coefficients of the gravitational potential 

terms is changed. To show that this may be true also for the angular component 

we list the results of calculations for an Earth whose potential includes the zeroth 

second and third zonal harmonics, and finally, the zeroth, second, third and 

fourth zonal harmonics. 

* 

Case III. 


" 2 3 *i 

~ tL 1 + ^ ^ (3 cos 2 0 - 1) +y |-!j (5 cos 3 6 - 3 cos 0)J 


(59) 


Here, 


d 2 = 

V r(/i + J, + j a > “ -■ 


mr g S 


zdt /i-i/i«y r 3 j l( 3 

1 sin 2 6 l L 


cos 2 6 - 1) + 20 J. 


© 


|cos 3 6 - ^ cos . 


(60) 


*(m + j 2 + j 3 ) 


^{' + }(if) pJ i <^o s =-.).20J a ^(co 5 3 9 -|co^)]J 


( 61 ) 
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Also, 


(f) [ 2J * * 4 (r)( 5 cos “e - C„! ej' 


(63) 


^<m + J 2 


$ [ 2 J» % * J . (r) ( 25 cos ‘ V 8 e t- - srj 


8' 

(64) 


and 


+ J 3 ) = COS 6 ^ 


2J_ + J 3 f-)f 25 cos 6 P ^3 \ — ^ 

2 3 \ r /\ fc cos cos 3 ft V 


2J ** J »^)( 5co * e s-shg) 


1/2 




>2 ’ J3 2 
(65) 


Case IV. 


= 7 [ l *r(r) (y) (■ 

- ^ 

+ i J 4 f— ' \ (35 cos 4 0 - 30 cos 2 0 + 3)1 


cos 3 6 - - cos 


(60 





/ 


♦SO.M..J,. i,> - -■£ (r) (t) ( s cos *g - 




(140 cos 2 <9,, - 60) 


3/x cos 2 6 / r \J 


j2+J 3 + ^> - 2rar 3- V7/ ‘ J2 + J 3 Vr/i — e "^7^" — T7 


l _ ) pj 2 + J3 V - /M® cos “ 




(140 cos 2 0 fe , - 60) . 


^(m + J 2 *j 3 + j 4 ) 55 cos 


2 J 2 + J3 (tJI 25 cos 0 e ~ 


£ cos 3 0g 


2J2 + J 3 u)( 5 COS ^fc 1 “ cos7 


+ ^ (140 cos 2 0g - 60) 


m 


(140 cos 2 e - 60) 


^+J 2 + J 3 + J 4 ) 


From the general equations (35) and (40, one can obtain £ r /u> x and <£ 0 /x 0 for 
any number of zonal harmonics. It is of interest to compare the coefficients of 
or the brackett terms of <j> 2 and a; 2 with the expression for the anomalistic 


mean motion, 
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where n 0 = (/x/a 3 ) x/7 for central force motion, a is the semi-major axis, and e 
is the eccentricity. The right hand side is expanded through the coefficient of 
the second zonal harmonic only, and can be compared to equations (50) and (51) 
for <b ? and*; 2 . 


V. TIME AVERAGES 

Figures 6(a) and (b) are pictorial representations of translational mechanical 
and series electrical circuits. The general differential equations for these 
systems are 


— — K. + ^ — + ky = F a cos cot, translational mechanical. 
a +2 dt 0 


(b) Ll f + R — + 1 = - o;E n sin cot, scries electrical. 
dt 2 dt c 


The vertical displacement coordinate of the weight is y, M is the vibicting 
mass, 'i is the friction, k, the spring constant, F 0 and ^ F , the amplitude and 
frequency of the impressed driving force, respectively. For the electrical 
case, L is the inductance, c, the capacitance, R, the resistance, E 0 and the 
amplitude and frequency of the impressed voltage, and i, the current coordinate 
flowing around the loop. If ^ (R) is zero, then the system is a conservative one. 
If no driving forces are present, then the equations are homogeneous, and the 
natural or resonant frequencies are given by 
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i/E. 

f M 


and ■■ - 

M vLC 


(75) 


The frequency be made to vary by changing the equivalent mass or spring con- 
stant of the system. In our problem, the mass is constant but the restoring force 
will depend upon the number of harmonics in the potential, or in other words, 
the distribution of the Earth's mass. From our results above, if we take J 2 to 
be approximately J 2 = -1.0825 x 10~ 3 , J 3 - 2.5450 x 1CF 6 , and J 4 = 1.6715 x 10“ 6 
where ^ = 1, then the effect of including the third and fourth zonal harmonics 
is to shift the value of o> r and oj $t opposite to thut of J 2 . However, since J 2 is of 
the order of 10 3 times J 3 or J 4 , the latter effect is very small. In any event, 
equations (35) and (40) predict a phase shift due to uncertainties and variations 
in the terms of the gravitational potential. The amplitude of the oscillation for 
a free, conservative system, representing a solution of equations of the type (74) 
is given by 

N = Njj cos(^t + b 0 ) t (76) 


for the cross-track error, and 


R = Ro cos t + 8 f ) % 


(77) 


for the radial uncertainty. Here N 0 and R 0 are the maximum error amplitudes 
for each case, and & is a corresponding phase a.'gle. Since the epoch times are 
taken to be the same for both observed and computed d^ta, then the phase angle 
is set equal to zero. Let us now examine the time averages of the ^ross and 
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radial uncertainties. The time averages of equations (76) and (77) are obtained 
by integrating both sides with respect to t from 0 to r , and dividing by r : 




cos t d t 


1 

T 


No 

°>e 


[cos (oj 6 t ) 1] 


(78) 


Since N Q and u> 9 are constant, when u Q r - 2n77, where n = 0, 1, 2, . . . , then, 


1_ 

T 



Ndt = 0. 


(79) 


That is, if r is chosen to be the period, then (78) vanishes. Even if uu q t is not 
an integer multiple of 2 77, the maximum value that the bracket portion of t?8) can 
be is 1, but choosing r sufficiently long, the right hand side of (78) can be made 
as small as desired. In this case, 


Similarly, 


fu-L. r 

T — 0 C J 
•>0 


Ndt -• 0 


/• T 

r = — - — r 

7- -*cn 1 

Jo 


Rdt - 0 


(80) 


(81) 
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VI. CONCLUSIONS 

In section in, we obtained relationships between the behavior of the radial and 
cross track errors and uncertainties in the orbital elements. In section IV, we 
found the equations relating this behavior to certain terms of the Earth's gravi- 
tational field. An examination of figures 2, 3, 4, and 5, indicates that the tracking 
errors oscillate harmonically with a frequency close to that of the orbital fre- 
quency. The wave interference analysis, shows that the radial frequency is one 
half the sum of the observed and computed orbital frequencies. These frequencies 
are functions of the semimajor axis of the orbit. The analysis of section IV shows 
that except for the case of a spherical Earth, the radial frequency is a function of 
the gravitational potential energy and the orbital frequency. Both results are 
consistent in that they show a proper frequency or phase shift for the radial 
oscillation. Results for the cross track case are similar, except that this 
phase shift does not appear until the computation potential includes the third 
zonal harmonic or higher. This indicates that quantitatively, there appears to 
be a lack of uncertainty in the eccentric anomaly and argument of perigee for 
the spherical Earth case, which is reasonable, since the orbit is then bound and 
closed and completely determined with full accuracy for all time. 

From these considerations, there emerges the following physical picture 
for the tracking error problem: The Earth-satellite system and the corresponding 
orbit dL Terences are considered as a harmonic oscillator system whose 
frequency is determined by the ratio of the spring constant to the mass for the 
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case of a translational mechanical system, or by the inverse product of the 
inductance (inertia) times the capacitance (elastance) for an electrical series 
circuit. In our case, the spring constants or capacitances are described in 
terms of the first and second spatial derivatives of ihe Earth's potential. The 
gravitational potential, which as a solution of Laplace's equation, involves here, 
only the ordinary Legendre polynomials, and these in turn, are determined by the 
distribution of mass in the Earth. This distribution therefore is described by the 
assigned values of the coefficients of these zonal harmonics. 

The results we have obtained in this paper arc valid for conservative systems, 
namely those in which the total energy and angular momentum are conserved. 

The frequency equations that we developed in section IV, assume an azimutlial!y 
symmetric potential. Because of this, the corresponding canonical momentum 
is a constant of the motion, and this led us to our analytic expressions for the 
tracking frequencies as functions of the potential. An important technique in 
theoretical physics is to look for symmetries within the system. For each 
symmetry, there follows a conservation law, which simplifies calculation of 
physical properties of the system. This approach, combined with the method 
of solution by separable Hamiltonian, makes it possible to do the gravitational 
theory cf a satellite orbit very accurately, and without the use of perturbation 
theory (References 5, 6). In addition, a method has been developed (Reference 7) 
which permits separation of the Hamilton Jacobi equation in the presence of a 
nonconservative force such as air drag. In this case, the atmospheric density 
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is incorporated into the equations of motion in an analytic fashion by fitting an 
exponential form to an accurate density profile. This allows one to integrate 
the variational equations of the orbital elements due to air drag in terms of 
elementary functions. This then has the effect of advancing the boundary condi- 
tions on the Hamilton-Jacobi equation. As a result, it is possible to include 
atmospheric resistance in orbital error analysis. 

It is expected that the results obtained here especially those of sections HI 
and IV, are to be extended in forthcoming analyses to more detailed cases such 
as given in figures 4 and 5, for GEOS-II. 
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FIGURE CAPTIONS 

Figure I* Two waves of nearly equal frequency (a) and with equal magnitudes, 

(b) are differenced, (c) to give a wave whose amplitude (dashed line) var- 
ies periodically. The effeec described from (b) to (c) is commonly 
referred to as the phenomenon of beats. 

Figure 2(a). Short term view of the radial error for the Cannonball satellite. 

Figure 2(b). Expanded time scale for the radial error of the Cannonball satellite. 

Figure 3. Expanded time scale of the cross track error for the Cannonball 
satellite. 

Figure 4. Tracking errors for the GEOS-II satellite using GEM 1 (doppler data). 

Figure £. Tracking errors for the GEOS-II satellite using 69 SAO Standard Earth 
(doppler data). 

Figure 6. Physical systems whose natural frequencies are equivalent to tracking 
error oscillations. 







Figure 1. Two waves of nearly equal frequency (a) and with equal magnitudes, (b) are dif- 
ferenced, (c) to give a wave whose amplitude (dashed lino) varies periodically. The 
effect described from (b) to (o) is commonly referred to as the phenomenon of beats. 
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Figure 3. Expanded time scale of the cross track error 
for the Cannonball satellite 
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Figure 6. Physical systems whose natural frequencies are equivalent 
to tracking error oscillations 
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